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^ ■ Abstract 

Q 

\^ • We consider the inverse scattering problem on the energy interval in three dimensions. 

PsJ , We are focused on stability and instability questions for this problem. In particular, 

we prove an exponential instability estimate which shows optimality of the logarithmic 
' stability result of [Stefanov, 1990] (up to the value of the exponent). 

< 

'5 : 1. Introdution 

We consider the Schrodinger equation 

- Alp + v{x)iIj = Eijj, xGM^ (1.1) 

where 



;i.2) 



> 

V/-^ . V is real- valued, v G L°°(R 

I v{x) = 0{\x\~^~'^), \x\ — > oo, for some e > 0. 

CN ■ 

Under conditions (11.21) . for any /c G \ equation (II. ip with E = k"^ has a unique 
O ■ continuous solution (x, k) with asymptotics of the form 

"O : r{x, k) = e^'^ - 27r''—-f A' 1^1 + 4 n 

|x| \ |/c| / \ \x\ 

S : (1-3) 

(X 
uniformly in - — - 
\x\ 

where f{k/\k\,u, \ k\) with fixed is a continuous function oi u E S"^. 

The function f{6,u,s) arising in (II. 3p is refered to as the scattering amplitude for the 
potential v for equation (II. ip . (For more information on direct scattering for equation (II. ip . 
under condition (11.20 . see, for example, [6j and [llj.) 

It is well known that for equation (II. ip . under conditions (II. 2p . the scattering amplitude 
/ in its high-energy limit uniquely determines -0 on , where 

v(p) = (27r)-^ / e'P''vix)dx, p G M^ (1.4) 

via the Born formula. As a mathematical theorem this result goes back to [5] (see, for 
example. Section 2.1 of |TT] and Theorem 1.1 of [13] for details). 
We consider the following inverse problem for equation (II. ip . 
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Problem 1.1 Given / on the energy interval J, find v. 

In [7j it was shown that for equation (II. ip . under the conditions ( II. 2p . for any E > and 
6 > the scattering amphtude f{6,u,s) on {{6,u,s) e S"^ x S"^ x E < s"^ < E + 6} 
uniquely determines v{p) on {p G | \p\ < 2\/E}. This determination is based on solving 
linear integral equations and on an analytic continuation. This result of [7] was improved 
in On the other hand, if v satisfies (II. 2p and, in addition, is compactly supported or 
exponentially decaying at infinity, then v{p) on {p e | |p| < 2\/E} uniquely determines 
v(p) on {p G I \p\ > 2\/E} by an analytic continuation and, therefore, uniquely determines 
V on M^. 

In the case of fixed energy and potential v, satisfying (II. 2p and, in addition, being 
compactly supported or exponentially decaying at infinity, global uniqueness theorems and 
precise reconstructions were given for the first time in [12j, [13j. 

An approximate but numerically efficient method for finding potential v from the scattering 
amplitude / in the case of fixed energy was devoloped in [15j. Related numerical implementation 
was given in [2j. 

Global stability estimates for Problem 1.1 were given by Stefanov in [17j (at fixed 
energy for compactly supported potentials), see Theorem 2.1 in Section 2 of the present 
paper. In [17], using a special norm for the scattering amplitude /, it was shown that the 
stability estimates for Problem 1.1 follow from the Alessandrini stability estimates of [Ij for 
the Gel'fand-Calderon inverse problem of finding potential v in bounded domain from the 
Direchlet-to-Neumann map. The Alessandrini stability estimates were recently improved by 
Novikov in [TBj . 

In the case of fixed energy, the Mandache results of [10] show that logarithmic stability 
estimates of Alessandrini of [1] and especially of Novikov of |T6J are optimal (up to the 
value of the exponent). In [H] studies of Mandache were extended to the case of Direchlet- 
to-Neumann map given on the energy intervals. Note also that Mandache-type instability 
estimates for the elliptic inverse problem concerning the determination of inclusions in a 
conductor by different kinds of boundary measurements and the inverse obstacle acoustic 
scattering problems were given in [3]. 

In the present work we apply to Problem 1.1 the approach of [10], [8] and show that the 
Stefanov logarithmic stability estimates of [17] are optimal (up to the value of the exponent). 
The Stefanov stability estimates and our instability result for Problem 1.1 are presented and 
discussed in Section 2. In Section 3 we prove some basic analytic properties of the scattering 
amplitude. Finally, in Section 5 we prove the main result, using a ball packing and covering 
by ball arguments. 

2. Stability and instability estimates 

In what follows we suppose 

snppv{x) C D = B{0,1), (2.1) 

where B{x,r) is the open ball of radius r centered at x. We consider the orthonormal basis 
of the spherical harmonics in L'^{S'^) = L^[dD): 

{Yj':j>0; l<p<2j + l}. (2.2) 

The notation {aj^p^j2P2) stands for a multiple sequence. We will drop the subscript 

< Ji, 1 < Pi < 2ji + 1, < j2, 1<P2< 2j2 + 1. (2.3) 
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We expand function f{6,uj, s) in the basis {Yj'_^^ x Yj'^^}: 

f{e,u,s)= ^nv^nvMyf^mSi^)- (2.4) 



JliPlJ2,P2 



As in 17 we use the norm 



2ji+2cri . , N 2j2+2cr2 



ll/(-,-,^)lk.2 = <| E (^) (^) l«..™2(^)P^ • (2.5) 



1/2 



If a function / is the scattering amphtude for some potential v G L°°{D) supported in 
B{0,p), where < p < 1, then 

N ji+3/2 / \ j2+3/2 

esp \ ' ( esp ^ ' 



l%mi2P2(s)l < C{s,\\v\\l^^d)) [2j^:p[ j \ 2j2 + l ) ^^'^^ 

and, therefore, ||/(-, ■, ■5)||o-i,o-2 < oo? see estimates of Proposition 2.2 of [17j. 

Theorem 2.1 (see [IZ]). Let Vi, V2 be real-valued potentials such that Vi G L°°{D) fl H'^(M.^), 
suppvi C -8(0, p), I < N for i = 1,2 and some N > 0, q > 3/2 and < p < 1. Let 

fl and f2 denote the scattering amplitudes for vi and V2, respectively, in the framework of 
equation U.l\) with E = s"^, s > 0, then 

\\Vl - V2\\l^(D) < c{N,p)(j)s{\\fi{-,-,s) - /2(-, -,5)113/2,-1/2), (2.7) 

where 4>s(t) = ( — Int)^"^ for some fixed 6, where, in particular, < S < 1, and for sufficiently 
small t > 0. 

The main result of the present work is the following theorem. 

Theorem 2.2. For the interval I = [si, S2], such that Si > 0, and for any m > 0, a > 2m 
and any real ai,a2 there are constants /3 > and N > 0, such that for any Vq G C"^{D) 
with ||fo||L°°(D) ^ N> suppvo C 5(0,1/2) and any e G (0,A^), there are real-valued potentials 
Vi,V2 G C"^{D), also supported in B{0, 1/2), such that 

sup(||/i(-,-,s) - /2(-, -,5)11^1,^2) < exp (-e"^) , 

\\v\ - V2\\l^(d) > e, (2.8) 
\\vi - vo\\l'^{d) < e, i = 1,2, 
\\vi - follc'Hi?) < /5, ^ = 1,2, 

where /i,/2 are the scattering amplitudes for Vi,V2, respectively, for equation ( fi.ij) . 

Remark 2.1. In the case of fixed energy Si = S2 we can replace the condition a > 2m in 
Theorem 12.21 by a > 5m/ 3. 

Remark 2.2. We can allow /3 to be arbitrarily small in Theorem 12.21 if we require e < eo 
and replace the right-hand side in the first inequality in (12. 8p by exp(— ce""), with eo > 
and c > depending on /3. 

Remark 2.3. Note that Theorem 12.21 and Remark 2.1 imply, in particular, that for any real 
(Ji and (J2 the estimate 

||fi -V2\\l^(d) < c(iV, p,m,/) sup 05(1 |/i(-, -,5) - /2(-,-,s)||^,.^J (2.9) 

sei 
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can not hold with 6 > 2m in the case of the scattering amphtude given on the energy interval 
and with 6 > 5m/ 3 in the case of fixed energy. Thus Theorem 12.21 and Remark 2.1 show 
optimality of the Stefanov logarithmic stability result (up to the value of the exponent). 
Remark 2.4. A disadvantage of estimate ( 12. 7p is that 

6 < 1 even if m is very great. (2.10) 

Apparently, proceeding from results of [16j, it is not difficult to improve estimate (12. 7p for 

6 = m + o{m) as m — 7- oo. (2-11) 

3. Some basic analytic properties of the scattering amplitude 

Consider the solution tp^{x, k) of equation 11.11 see formula (II. 3p . We have that 

V'+(x,A;) = e^^>+(x,^,s), (3.1) 
where 6 G S"^, k = s9 and fJ^^ix, 6, s) solves the equation 

^+(x, e,.s) = l- [ G^ix, y, s)e-'^'^^-y^viy)fi^iy, 6, s)dy, (3.2) 



where 



^is\x-y\ 



G+(x,i/,s) = — (3.3) 

4:'k\x — y\ 

We suppose that condition (12. ip holds and, in addition, for some /i > we have that 

|Ims|</i, (3.4) 



where D = B(0, 1) 



Then, in particular. 



ci{h,D)\\v\\Lo^^D)<l/2, (3.5) 

Ci(/i, D) = sup / — j -dy. (3.6) 

xeD Jd 4vr|x - y\ 



^^-ise{x-y) ^is\x-y\\^ < g2h|a;-j/| -j?^ 

Solving ( 13. 2p by the method of succesive approximations in L°°{D), we obtain that 

\fi+{x,9,s)\<- ^ , eeS^, xeD. (3.8) 

Lemma 3.1. Let cijipij2P2{^) denote coefficients f{s, 6, u) in the basis of the spherical harmonics 
{yj^^ ^^^^}' "^here f is the scattering amplitude for potential v G L'^{D) such that conditions 
and 5\) hold for some h > 0, 

f{e,u;,s)= <^nP.npMyf:mj:\u;). (3.9) 

jl,Pl,j2,P2 

Then cij^p^j2P2{s) is holomorphic function in Wh = {s \ \Ims\ < h} and 

l«iiPU2P2(5)l < C2{KD)s G Wh. (3.10) 
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f{e,uj,s) = I e''^^-^^''v{x)n+{x,e,s)dx. (3.11) 



Proof of Lemma \3.1\ We start with the well-known formula 

1 

Note that, since 6,uj G S*^, 

|gis(6»-w)a;| ^ ^2\Ims\\x\ j^g ]^2) 

Combining it with ([2J]), ([33]), ([33]) and dXTT]) we obtain that 

\f{e,uj,s)\ < C2ih,D) for s G (3.13) 

Using also that 

anP.npM= [ fie,u,s)Y[;{e)Yf^'iu)dedu (3.14) 
we obtain the result of Lemma 13.11 ■ 



A fat metric space and a thin metric space 



Definition 4.1. Let (X, dist) be a metric space and e > 0. We say that a set y C X is an 
e-net for Xi C X if for any a; G Xi there is ?/ G F such that dist{x, y) < e. We call e-entropy 
of the set Xi the number 'He(Xi) := log2min{|y| : F is an e-net fot Xi}. 

A set Z C X is called e-discrete if for any distinct zi, Z2 G Z, we have dist{zi, Z2) > e. We 
call e-capacity of the set Xi the number := log2max{|Z| : Z C Xi and Z is e-discrete}. 

The use of e-entropy and e-capacity to derive properties of mappings between metric 
spaces goes back to Vitushkin and Kolmogorov (see [9j and references therein). One notable 
application was Hilbert's 13th problem (about representing a function of several variables 
as a composition of functions of a smaller number of variables). In essence. Lemma [4.11 and 
Lemma [4.21 are parts of the Theorem XIV and the Theorem XVII in [9]. 

Lemma 4.1. Let d > 2 u m > 0. For e, /3 > 0, consider the real metric space 

XmeP = {V e C"^{R'^) I SUppV C 5(0, 1/2), | |t; 1 1 Loomd) < e, \\v\\c^r^d) < f3} 



with the metric induced by L°°. Then there is fi > such that for any /3 > and e G (0, fi(3), 
there is an e-discrete set Z C X^e^ with at least exp (^2^'^^^(fi(3 / e)'^^"^^ elements. 

Lemma 4.2. For the interval I = [a, b] and 7 > consider the ellipse Wi^^ G C.- 

,a + b a — b , , ^ , . , . 

^/,7 = {—^ \ ^cos2; I \Imz\ < 7}. (4.1) 

Then there is a constant v = u^C, 7) > such that for any 6 G (0, e~^) there is a S-net for 
the space of functions on I with L°°-norm, having holomorphic continuation to Wj^^ with 
module bounded above on Wj^^ by the constant C, with at most exp(i/(ln 5~^)^) elements. 
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Remark 4.1. In the case of a = 6, taking 



(4.2) 



we get 5-net with at most exp(i^ln5~^) elements. 

Lemma [4. II and Lemma [4.21 were also formulated and proved in |10) and [8], respectively. 

For the interval / = [si,S2] such that Si > and real (Ji,(J2 we introduce the Banach 
space 

(4.3) 



X 



< oo 



where 

%lPli2P2('^y 



-^1,^1, 



sup 



2ji + lV'^"' /2j2+ 1^''+"' 



es 



es 



{s)\\. (4.4) 



We consider the scattering amplitude / for some potential v G L°°{D) supported in 5(0, p), 
where < p < 1. We identify in the sequel the scattering amplitude f{s, 6, u) with its matrix 



''3lPl32P2\ 



s) j in the basis of the spherical harmonics {Y^^ x Y^^} . We have that 



SUp||/(-,-,s)||^,,^2 < C3 



^jlPlj2P2\^ , 



(4.5) 



where (Xi — a2 = cr2 — = 3 and C3 = c^{I) > 1. We obtain (14. 5 p from definitions (12. 5p . (14. 4 p 
and by taking C3 > 1 in a such a way that 



E 



Jl,Pl J2,P2 



2ji + lVV2j2 + l' ' 



es 



es 



< C3. 



(4.6) 



For h > Owe denote by Ah the set of the matrices, corresponding to the scattering amplitudes 
for the potentials v G L°°{D) supported in -8(0, 1/2) such that condition (13. 5p holds. 



Lemma 4.3. For any h > and any real ai, a2, the set Ah belongs to Xj^, 



In addition, 



there is a constant rj = ri{I, h, o"i, (72) > such that for any 5 G (0, e ^) there is a 6-net Y 
for Ah in Xi^o-i,a2 with at most exp {rj (ln5~^)^ (1 + lnln(5~^)^j elements. 



Proof of Lemma \4-^ We can suppose that ai,(T2 > as the assertion is stronger in this 
case. If a function / is the scattering amplitude for some potential v G L°°{D) supported in 
5(0, 1/2), we have from ([2l]) that 



2ji + 1 



2j2 + 1 



i2+0"2 



I %lPli2P2 



»l < C4 



es / \ es 
where C4 — C4 (/, h) > 0. Hence, for any positive ai and a2, 



2^+32 



(4.7) 



%lPlj2P2 ('^) 



and so the first assertion of the Lemma 14.31 is proved. 



(2ji + l)-^(2j2 + l)-^ . ^ 

< sup 1 C4 I < 00 



(4.8) 
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Let /<5,cri,(T2 be the smallest natural number such that 04(2/ + ' < 5 for any 

^ > h,ai,a2- Taking natural logarithm we have that 

- Inc4 - (ai + (T2) ln(2/ + 1) + /ln2 > In^^ for / > /vi,^2- (4-9) 

Using ln6~^ > 1, we get that 

Vi,-2<C'lnr\ (4.10) 

where the constant C depends only on h, ai, o"2 and / = [si,S2]. We take Wj = Wj^^ of 
(14. ip . where the constant 7 > is such that Wj C {s \ \lms\ < h}. If max(ji,j2) < h,ai,a2y 
then we denote by Yj-^p-^^j^p^ some Sj-^p-^j^p^-net from Lemma 14.21 with the constant C = C2, 
where the constant C2 is from Lemma [3.11 and 



Otherwise we take l^ipij2P2 — {0}- We set 

^ ~ I (^'^iiPii2P2("^)^ I ^jiPij2P2{^) ^ ^jipij2P2^ ■ (4-12) 
For any ^ajipij2P2("5)^ ^ -4.^ there is an element (^&jipij2P2('5)) ^ ^ such that 
/2j^y^+- /2j2+lV^+'^^ I (\-h Ml< 

I I I eS / \^jrpij2P2\^ ) "jlPlj2P2\^)\ — 



(4.13) 



in the case of max(ji, J2) ^ h,ai,a2 ^'^^ 

'2ji + l\ /2j2 + l^''+"' 



I \ eS / \^hPlj2P2{^) ^jlPlj2P2{^)\ — 

{2j, + ir-{2j, + ir^ ^ (2max(ji,j2) + l)"^''"^ , , 
< C4 — — < Cd r- — ^1 < 0, 



(4.14) 



otherwise. 

It remains to count the elements of Y. We recall that |^ipij2P2l = 1 in the case of 
max(ji,j2) > h,cn,a2- Using again the fact that ln5~^ > 1 and ( I4.1QP we get in the case of 
max(ji, 72) < I5 

\YnP.nP.\ < exp(^(ln5^i,,,,J^) < exp (z/' (inr^)' (l + Inlnr^)') . (4.15) 

We have that ns^ai,a2 < ^5,0-1,^2 (2^'5,<7i,<72 + 1)^ < (2^5,^1,^2 + 1)^, where ns,aua2 is the number of 
four-tuples (ji,Pi, J2;P2) with max(ji, J2) < h,ai,a2- Taking rj to be big enough we get that 

|r| < (exp(z.' (lnri)'(l + lnlnri)'))"'" " 

< exp (u' {In6-^Y (l + Inlnr^)^ (1 + 2C"lnri)^) (4.16) 

< exp (^7] {\n6-^Y (1 +lnlnr^)' 



Remark 4.2. In the case of Si = S2, taking into account Remark 4.1 and using it in fl4.15p 
and fl4.16p . we get 5-net Y with at most exp (^rj (ln5~^)^ (1 + In In 5"^)^ elements. 
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5. Proof of Theorem 2.2 

We take N such that condition (13. Sp holds for any ||f||L°°(_D) < 2A^ for some h > 0. By 
Lemma [4. 11 the set vq + Xmeis has an e-discrete subset vq + Z. Since e G (0, A^) we have that 
the set Y constructed in Lemma 14.31 is also (5-net for the set of the matrices, corresponding 
to the scattering amplitudes for the potentials v & vq + Xmei3- We take 6 such that 20^6 = 

exp ^— e^'^^, see (14. 5p . Note that inequalities of (12. Sp follow from 

\vo + Z\>\Y\, (5.1) 

where the set Y is constructed in Lemma [4.31 with ai = o"i + 3 and a2 = <^2 + 3. In fact, if 
\vq + Z\ > \Y\, then there are two potentials Vi,V2 Ev^ + Z with the matrices {^jxpij2P2{s)^ 

and (^&jipij2P2('^)) ' corresponding to the scattering amplitudes for them, being in the same 
-^/,(Ti,cr2-ball radius 5 centered at a point of Y . Hence, using (14. 5 P we get that 



sup||/i(-,-,s) - /2(-,-,s)||^,,„2 < C3 (^ajipij2p2(s)) - (h 

< 2C36 = exp (^~^ °) ■ 



'JlPlj2P2 [S, 



< 



(5.2) 



It remains to find /3 such that (15. ip is fullfiled. By Lemma H73l for some rja = Tia{I, ai, (72, a) > 




\Y\<exp(^7j (\n{2c3) + e"-)^ (^1 + In (ln{2c3) + e"^))^^ < exp (^ri^t 



3 



) . (5.3) 
Now we take 

/3 > 1^-^ max {N, r^'^'^^"^) . (5.4) 



This fulfils requirement e < in Lemma [4.11 which gives 

m 



+ Z\ = \Z\ > exp (2-^(/i/3/e)3/"^) > 

> exp (2-^(77;^/322m/^^)3/m^ > 



(5.5) 



This completes the proof of Theorem 12.21 

In the case of fixed energy Si = S2, using Remark 4.2 in ( 15. 3p . we can replace the condition 
a > 2m in Theorem 12.21 by a > 5m/3. 
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